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Abstract 

Ferroelectric domain switching in c-axis-oriented epitaxial Pb(Zro. 2 Tio. 8)03 thin films was studied 
using biased scanning probe microscopy tips. While linear and logarithmic dependence of domain 
size on tip bias and writing time, respectively, are well known, we report an additional linear depen¬ 
dence on relative humidity in the 28-65% range. We map out the switched domain size as a function 
of both the tip bias and the applied pulse time and describe a growth-limited regime for very short 
pulses and a nucleation-limited regime for very low tip bias. Using “interrupted-switching” mea¬ 
surements, we probe the nucleation regime with subcritical pulses and identify a surprisingly long 
relaxation time on the order of 100 ms, which we relate to ionic redistribution both on the surface 
and within the thin film itself. 

Keywords 

ferroelectrics, polarization switching dynamics, piezoresponse force microscopy, humidity, nucle¬ 
ation 

1 Introduction 

The characteristic electric polarization of ferroelectric materials can be locally controlled via an 
electric field applied with a biased scanning probe microscopy (SPM) tip [1,2]. The extremely small 
domains - distinct regions with a well-defined polarization state - formed in this fashion remain sta¬ 
ble over extended time periods [3]. Such nanoscale non-volatile domains allow ultrahigh (beyond 
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Figure 1 . (a) Schematic of the electric field induced nucleation and subsequent growth of a domain in a ferroelectric thin 
film. Representative images of (b) topography, (c) PFM phase, and (d) PFM amplitude signals on our Pb(Zro. 2 Tio. 8)03 
(PZT) sample. 


10 TBit inch“^) information storage densities to be envisaged [4], motivating extensive studies of 
their switching dynamics [5]. Thermodynamically, the switching process can be described [6,7] as 
the initial, almost instantaneous, softening and reorganization of a highly localized volume of polar¬ 
ization under the biased SPM tip, activated by the intense electric field in this region, leading to the 
formation of a critical nucleus with charged domain walls, then very rapid forward propagation, and 
finally slower lateral expansion of the growing domain, as schematically illustrated in figure 1(a). 
Experiments have focused primarily on the latter, more accessible regimes, using piezoresponse 
force microscopy (PFM) to image needle-like [8] or cylindrical [9] domains in single crystals and 
thin films. Such studies have demonstrated the creep-like lateral motion of domain walls pinned 
by defects [10-12], recently confirmed with atomic resolution by direct imaging of domain growth 
using transmission electron microscopy [13]. They also highlighted the crucial role of surface ad¬ 
sorbates [14-17], particularly water, which determine the electrostatic boundary conditions and thus 
the size, shape and growth rate of the domain [18-20]. In contrast, the subcritical regime, where the 
voltage pulses are too short to fully stabilize domains of inverted polarization, is much less well un¬ 
derstood, with the additional complexity of highly dynamic electrochemical effects such as charge 
injection or re-ordering under the biased SPM tip [21]. These different processes all take place in 
a very confined region of high electric field, at possibly very different time scales, and can strongly 
influence the initial stages of polarization softening and stabilization of the critical nucleus. 

In this letter, we present a detailed PFM mapping of domain switching in epitaxially-grown 
Pb(Zro. 2 Tio. 8)03 (PZT) as a function of tip bias, writing time, and relative humidity in the range 
from 28% to 65%. For low tip bias, we find that successful switching occurs only when the activa¬ 
tion threshold for polarization reversal is itself comparably low, and thus depends very sensitively 
on the local, highly heterogeneous defect distribution [22,23] which provides a randomly varying 
pinning/activation energy landscape [12]. However, when these activation threshold requirements 
are met, and as long as the tip bias is maintained, the resulting domains nonetheless grow outward 
by radial domain wall motion, up to significant size. For short writing times and high tip bias, in 
contrast, we observe a very high probability of initial nucleation, but successful switching depends 
on the stabilization of this critical nucleus under the SPM tip. The resulting domains are very small. 
We refer to these two regimes as nucleation-limited and growth-limited, respectively, in analogy to 
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switching behavior observed in ferroeleetrie eapaeitors [24] and thin films [25]. Finally, using very 
short repeated pulses well below switehing limits, separated by a defined “interruption time”, we 
probe the suberitieal regime of domain switehing, demonstrating that the initial polarization soft¬ 
ening and suberitieal nueleus formation are eharaeterized by surprisingly long relaxation times of 
over 100 ms. 

We used a 270 nm thiek epitaxial, c-axis oriented ferroeleetrie PZT thin film grown by olf-axis 
radio-frequeney magnetron sputtering on 35 nm thiek eondueting SrRuOs on (001) single-erystal 
SrTiOa (CrystTec), with high erystalline quality [26] and 2.4 nm root-mean-square surfaee rough¬ 
ness. The remanent polarization of 75 [xC om“^ is perpendieular to the film plane, and monodomain 
(‘up-polarized’) as grown, showing an imprint of -1.6 V with eoereive voltages of -1-3.3 V and 
-6.5 V at a eyele frequeney of 500 Hz. To switeh the polarization, voltage pulses were applied 
to Bruker MESP tips with the SrRu 03 layer grounded, using a Veeco Dimension V atomie foree 
mieroseope [27]. PFM imaging with simultaneous reeording of topography (figure 1(b)), phase 
(figure 1(e)) and amplitude (figure 1(d)) signals was performed in eontaet resonanee using a drive 
frequeney 10 kHz above the resonanee peak (~300 kHz) at a drive amplitude of 2.0 V. The elfeetive 
domain radius was extraeted from the total area obtained from the binarized PFM phase images. 
A home-made easing around the mieroseope, eombined with either desieeant material (siliea gel 
beads) or evaporative water souree allowed the relative humidity (RH) of the sample environment 
to be stabilized in the range between 28% and 65%. The ambient relative humidity of the laboratory 
was measured to vary between 30% and 54% during the eourse of one year. 

2 Results and Discussion 

To map out the switehed domain radius as a funetion of the writing parameters, we applied a tip 
bias between 2 and 10 V for writing times ranging from 80 qs to 600 s. For eaeh set of parameters 
measured, represented by a blaek dot in figure 2, five separate switehing attempts were made, and 
the size of the resulting sueeessfully switehed domains averaged to obtain the elfeetive radius. Sue- 
eessful switehing was defined as a 180° PFM phase eontrast between the two polarization states, 
and the presenee of a well-defined ring feature eorresponding to a deerease of PFM amplitude at 
the domain wall, eharaeteristie of a eolumnar domain penetrating through the full film thiekness. 
At both RH values, we reeognize the previously-reported [12] linear (figure 2(e)) and logarithmie 
(figure 2(d)) dependenee of domain size on tip bias and writing time, respeetively. However, for 
eomparable ehoiee of writing parameters, domain radii at 44% RH (figure 2(b)) are almost twiee 
as big as those obtained at 29% RH (figure 2(a)). The minimum stable domain radius, almost in¬ 
dependent of the tip bias for an appropriate ehoiee of pulse time, is around 17 nm for 44% RH 
and only 10 nm for 29% RH. These very small domains remained stable for at least 6 weeks. The 
hatehed area represents the writing parameters for whieh the sueeess rate of domain switehing is 
less than 100%. In this region of parameter spaee, eomparable for both high and low RH, there is a 
progressive transition from eertainty of sueeessful switehing, with the resulting domain penetrating 
through the film thiekness, to a regime where the switehing and stabilization of a domain beeome 
rare. We note that unsueeessful switehing ineludes observations of no ehange in either phase or 
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Figure 2. Average effective domain radius as a function of tip bias and pulse time obtained at a relative humidity of 
(a) 29 % and (b) 44%. Iso-size lines and color shading indicate the parameters for which domains have the same radius. 
Hatching represents switching success rates below 100%. No switching was observed in the white region. Black 
dots represent each measured set of parameters, (c) Cross-sections along the white horizontal dotted lines through the 
parameter space showing the linear dependence of domain radius on tip bias, (d) Cross-sections along the white vertical 
dashed lines through the parameter space showing the logarithmic dependence of domain radius on pulse time. 


amplitude in subsequent PFM imaging, and of traee features, sueh as ean be seen in the upper row 
of domains in figure l(e, d), whieh nonetheless fail to meet the full eriteria for sueeessful switehing. 

For high tip bias but short writing time, the latter appears to be the key limiting faetor. In 
this region of parameter spaee, partially switehed needle-like domains [28], with eharged domain 
walls, whieh do not penetrate through the full sample thiekness ean often be deteeted, disappearing 
between hours and a few days after writing. Sueh partially switehed domains are eharaeterized by 
an intermediate PFM phase eontrast of less than 180° and a lower PFM amplitude signal without 
a distinet ring feature. Their presenee demonstrates that here stabilization and subsequent growth 
of the eritieal nueleus, rather than the initial nueleation step, limit domain switehing. For low tip 
bias and long writing time, in eontrast, we deteet no traee of partial switehing. Rather, any domains 
observed are fully stabilized, sometimes with diameters well above 100 nm, although the probability 
of their formation gradually approaehes 0%. The limiting faetor in this region of parameter spaee is 
thus the low tip bias, insuffieient to meet the threshold requirement for nueleation. This nueleation 
threshold, shown to loeally vary by over 6 V in epitaxial PZT [22], is determined by the distribution 
of defeets in the thin film. The domain switehing probability thus depends very sensitively on the 
loeal defeet landseape. We note that the upper limit of the hatehed area does not follow the iso-size 
lines of equal domain radius, thus for a tip bias below 4.5 V, domain switehing is not eertain even 
at 100 s pulse times. 
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Figure 3. (a) Linear relation between the domain radius and the relative humidity for pulses of 0.1, 1, 10, or 100 s, at a 
tip bias of 10 V. (b) Domain radius as a function of relative humidity for pulses with a total time of 1 s at a tip bias of 
10 V. Pulses were applied singly, or split into 10, 100, 1000, or 10000 intervals with an interruption time of 60 ms. 


To better understand the specific contributions of the different parameters, including humidity, 
we applied pulses of 0.1, 1, 10, or 100 s with a tip bias of 10 V to the sample, and measured the 
effective radius of over 1500 domains across the available RH range [29]. As shown in figure 3(a), 
we find a linear increase of domain size with the relative humidity in the 28-65% range, with domain 
radii at high RH being more than twice as big as those at low humidity. For different pulse times the 
slope of this linear dependence varies from 1.6 to 5.3 nm %RH“^. Since the vast majority of PFM 
measurements are carried out at ambient conditions, this very significant effect should be recognized 
in any quantitative study of domain switching and growth, especially given the annual RH variation 
in a standard laboratory. We note that a simple extrapolation of the observed linear dependence of 
domain size on humidity would yield zero-size switching limits in the 5-8% RH range. However, 
domain switching has been successfully demonstrated even under ultrahigh vacuum [20,30]. We 
therefore expect a deviation from the linear behavior for lower humidities approaching 0% RH. For 
long writing times at high humidity, we also observe an increasing spread in the measured domain 
size. Radii of domains written with 100 s pulses below 30% RH present a standard deviation of 
only 3-10%, while for domains written above 60% RH a much greater standard deviation of up 
to 25% is extracted. This increasing spread suggests a greater influence of local variations in the 
defect landscape, which would be expected as the electric field intensity and thus domain wall speed 
decreases further away from the tip. 

A straightforward explanation of the influence of the relative humidity on the growth of ferro¬ 
electric domains is readily provided by the presence of a water meniscus at the point of contact 
between the SPM tip and the sample, as imaged by Weeks et al. [31]. However, the several-100- 
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nm-wide menisci reported by these authors were formed only after RH saturation, then progressive 
decrease to a lower target value. As long as the system stayed in a RH range below 60%, corre¬ 
sponding to standard ambient conditions, the size of the meniscus at the SPM tip in fact remained 
below the 50 nm imaging resolution of the experiment. Thus, a theoretical estimate of the menis¬ 
cus size obtained using the Kelvin equation for capillary condensation [32] provides a more useful 
approach. In our case, with the assumption of a 1 nm thick adsorbed layer of surface water on top 
of the sample [33], we calculate the horizontal extent of the meniscus to range from 12 to 16 nm 
for RH between 28% and 65% [34]. In electrostatic simulations using finite element analysis with 
COMSOL Multiphysics its presence indeed both increases the lateral extent over which high values 
of the vertical electric field (along the out-of-plane polarization axis) are maintained, and decreases 
the maximum field intensity at the tip-sample contact point. A water meniscus would thus lower 
the probability of nucleation, which depends on the intensity of the electric field in the local region 
under the biased SPM tip, but would promote more rapid domain growth through the high field 
region were nucleation to occur, in qualitative agreement with the experimental results. However, 
the calculated effect remains very local, with no significant change in the field profile at distances 
beyond 35 nm from the contact point [34]. In contrast, experimentally we observe that domain 
growth is still linearly influenced by humidity up to distances corresponding to the largest measured 
domain radii of 350 nm, suggesting that the meniscus size is only a part of the physical mechanism 
and that longer range effects need to be considered. 

One possibility proposed by Brugere et al. [35] is that the electric field itself slowly propagates 
through the weakly conductive surface water layer, driving lateral domain growth. However, Guy- 
onnet et al. [20] demonstrated qualitatively similar domain wall creep both in ambient conditions 
and under ultrahigh vacuum, disfavoring this hypothesis. Their work highlighted instead the im¬ 
portant screening effects of surface water with respect to both the ferroelectric dipolar forces, and 
to defects which can act as nucleation and pinning sites. Varying the strength of this screening in 
a simple Ginzburg-Landau switching model with numerically introduced disorder reproduced the 
order-of-magnitude difference in the size of the resulting domain with identical tip bias and writing 
voltage. We also note that under the locally very intense electric field of the biased SPM tip, the 
disorder potential itself can dynamically evolve. Electrochemical effects such as oxygen vacancy 
injection or (re)ordering, surface charging, and even irreversible physical damage to the sample [21] 
can significantly change the defect landscape during switching, particularly at ambient conditions. 
The resulting complex interplay between surface water screening and defect dynamics would influ¬ 
ence not only the lateral growth of domains, but also the initial stabilization of the critical nucleus. 

To further explore these effects we used “interrupted-switching” PFM, measuring domain 
growth after the application of 10 V bias to the SPM tip for a total time of 1 s, subdivided into 10, 
100, 1000, or 10000 shorter pulses separated from each other by an interruption time of 60 ms, 
during which the tip was grounded. Yang et al. [36] reported that a similar use of pulse splitting into 
shorter intervals had no effects on the growth of the resulting domain. However, as the constituent 
pulses approach the limit for switching, on the order of 10-100 ps in the case of our sample, we 
would in fact expect significantly reduced domain sizes and lower probability of successful switch¬ 
ing. Such short pulses should thus allow us to probe the initial stabilization of the critical nucleus 
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under different RH. 

As ean be seen in figure 3(b), subdivision of the total pulse time into intervals longer than 10 ms 
indeed has relatively little effeet. While the resulting domain size may be slightly smaller than that 
written with a single pulse, the linear dependenee of domain size on RH is maintained aeross the 
measurement range. However, we find that the role played by RH becomes less important when 
the total pulse time is subdivided into intervals shorter than 1 ms. In fact, for 10000 constituent 
pulses of 100 [.IS the final domain size of 65-70 nm appears to be independent of the humidity. 
We confirmed these observations in repeated experiments for total pulse times of 10 and 100 s, 
and also investigated the effect of varying interruption times (60, 120, 280, and 560 ms) which 
produced no significant change in final domain radius [37]. In all these measurements with 10 V 
tip bias, the probability of successful switching was almost 100%, even for the shortest subdivision 
times. In contrast, we find that lowering the tip bias very strongly decreases this probability. In 
identical interrupted-switching PFM measurements with 8 and 6 V tip bias, we again observe that 
the subdivision of the total pulse time into many very short intervals results in decreased domain 
size, with little or no RH dependence. However, the probability of successful domain switching 
after 10000 pulses of 100 [xs is only 45% at 8 V, and the value drops to 15% for 6 V. For the lower 
(6 V) tip bias, domain writing using fractional pulses is also much more sensitive to tip wear and 
higher humidity levels [37]. 

These results, in line with the data presented in figure 2, confirm that nucleation and full sta¬ 
bilization of a narrow columnar domain under the biased SPM tip occur at relatively short time 
scales on the order of 100 [xs. However, for the subsequent lateral growth and stabilization of large 
domains, we find that both high RH and longer continuous exposure to the electric field are needed. 
In addition, for similar pulse times and lower tip bias the resulting written domains are not only 
smaller, but also characterized by decreased switching probability, pointing to the strong influence 
of local variations of the defect landscape. 

Finally, to access the subcritical regime of the switching dynamics we used interrupted- 
switching PFM with constituent pulses below the switching limit, with varying interruption times. 
Each individual 10 V pulse was only 10 [xs long, 10 times shorter than the threshold pulse duration 
necessary to successfully switch a domain. As shown in figure 4, in each case we wrote a domain 
array, doubling the number of single pulses at each line, starting at 1 (top) and ending at 512 
(bottom). The interruption times ranged from 60 ms (figure 4(a)) to 330 ms (figure 4(d)). For single 
10 [.xs, 10 V pulses we detect no visible effect on the ferroelectric film during subsequent imaging 
in either PFM phase or amplitude. However, with a repeated application of as few as 16 such 
pulses separated by interruption times of up to 100 ms, successful domain switching can be ob¬ 
served, presenting the minimum PFM amplitude ring feature characteristic of well-defined domain 
walls [38]. Moreover, even down to 4 repeated pulses, we can observe traces of partial switching, 
with intermediate PFM phase contrast and a somewhat amorphous PFM amplitude drop at the 
position where the bias pulses were applied. As presented in figures 4(e, f), longer interruption 
times decrease both the complete and partial switching success rates, determined by counting how 
many domains out of the 10 in each line were fully stabilized or partially switched. For interruption 
times of 170 ms, repeated application of 256 and 512 voltage pulses still allowed successful domain 
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Figure 4. PFM phase images of domains written with multiple pulses of 10 ps at 10 V tip bias at 33% RH. The number 
of individual pulses (indicated at left) ranges from 1 to 512, and the interruption time is either (a) 60 ms, (b) 100 ms, (c) 
170 ms, or (d) 330 ms. The leftmost column of each array (framed by dashed box) was written with single 10 ms pulses 
and serves as marker to locate the lines. Success rates at which (e) complete domain switching and (f) traces of partial 
switching can be seen in the phase signal as a function of the number of individual pulses, for the above mentioned 
interruption times. 


switching, but no stable domains were observed with 330 ms interruption times. Nonetheless, even 
with these very long interruption times, we still deteet a high rate of partial switching. 

The eumulative effect of individual pulses, which in themselves are too short to have any mea¬ 
surable effeet, proves that the as-grown state of the thin film is nonetheless modified after exposure 
to the eleetric field of the SPM tip for 10 p,s, rendering it sensitive to further exeitation. These highly 
loeal ehanges assoeiated with the formation and stabilization of the eritieal nueleus in the first stage 
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of switching persist for surprisingly extended periods. To fully stabilize a eolumnar domain through 
the film thiekness, the minimum total pulse time (160 p,s) is eomparable to the 100 p,s required with 
a single eontinuous pulse, suggesting that the ehanges indueed by eaeh eonstituent pulse are eon- 
served even with interruptions whieh are 10 000 time longer than the eonstituent pulses themselves. 
As the interruption time is inereased to 170 ms, however, the required total pulse time for the same 
result is 25 times longer, showing that in this ease the film returned at least partially towards its 
initial state during the interruptions. We thus estimate the relaxation time to be on the order of 
100 ms. 

This relaxation time is far beyond what would be expeeted for either eleetronie or even lattiee 
aeeommodation, or for rapid “baek-switehing” observed as domain walls rapidly adjust their posi¬ 
tion when the loeal eleetrie field of the tip is switehed off [13], and points to slower dynamie effeets. 
In addition, the formation and initial stabilization of the suberitieal nueleus appears to be largely RH 
independent within the measured range. One possibility is very loeal displaeements of and OH“ 
speeies in the immediate neighborhood of the tip [19,39], where the speeifie shape and extent of the 
water meniseus plays no signifieant role. These eharged speeies, formed at even the lowest partial 
pressures of water vapor on ferroeleetrie surfaees [40], have been previously shown to play a eritieal 
role in sereening and polarization switehing kineties [41]. Another key faetor is the redistribution 
of oxygen vaeaneies, for whieh drift veloeities of 5 pm s“^ under fields on the order of MV em“^ 
have been reported in PZT eapaeitors [42]. These defeets, eommon in all perovskite oxides, are 
partieularly abundant in thin films grown on SrTiOa substrates [43], and tend to reorganize under 
the influenee of a biased SPM tip [21]. Even small ehanges in the oxygen eontent in perovskite ma¬ 
terials have been reeently shown to have signifieant effeets [44,45] in terms of resistivity or e-axis 
parameter. A loealized inerease in the density of oxygen vaeaneies under the tip eould thus strongly 
eontribute to the observed long-life “aetivation” of this region and its subsequent tendeney to switeh 
more easily. 

We also note that reeent measurements of the role of RH in switehing on thinned LiNbOa single 
erystals, generally eontaining fewer defeets, showed eompletely opposite behavior with no observed 
ehange in domain size for 0-60% RH and a preeipitous drop thereafter [17]. These differenees 
strongly suggest that eleetroehemieal eontributions from within the ferroeleetrie, as well as ionie 
effeets of adsorbates on the surfaee, need to be eonsidered. 

3 Conclusions 

From a detailed mapping of the size of SPM-tip-written ferroeleetrie domains in PZT as a funetion 
of both tip bias and pulse time for whieh it is applied, we identify nueleation- and growth-limited 
switehing regimes at low bias and short writing times, respeetively. We show that these regimes 
depend on very different faetors, the first eontrolled by variations in the defeet landseape whieh 
loeally modify the nueleation bias threshold, the seeond by the dynamies of the stabilization of 
the eritieal nueleus and the initial stages of vertieal domain growth through the film thiekness. 
In addition, we quantify the linear dependenee of domain size on the relative humidity over and 
beyond the full range of values eorresponding to its annual variation in a standard laboratory, whieh 


9 


can significantly affect the measurements of domain switching dynamics in ambient eonditions. 
Finally, introducing interrupted-switching PFM in the suberitieal regime, we extraet surprisingly 
long relaxation times on the order of 100 ms for the loeal modifications in the ferroeleetrie film 
assoeiated with the formation and stabilization of the eritieal nucleus. 

Supporting Information 

Effeets of tip wear, geometry of the water meniseus, finite element simulations of the vertieal 
eleetrie field, alternative interrupted-switching PFM protocols, and PFM amplitude images of the 
suberitieal switching experiments. This material is available free of charge via the Internet at 
staeks. iop. org/njp/17/013002/mmedia. 
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